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A recent work showing that homogeneous and isotropic cosmologies involving scalar fields are 
equivalent to the geodesies of certain eflFective manifolds is generalized to the non-minimally coupled 
and anisotropic cases. As the Maupertuis-Jacobi principle in classical mechanics, such result permits 
us to infer some dynamical properties of cosmological models from the geometry of the associated 
effective manifolds, allowing us to go a step further in the study of cosmological dynamics. By 
means of some explicit examples, we show how the geometrical analysis can simplify considerably 
the dynamical analysis of cosmological models. 
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I. INTRODUCTION 



- - - PACS numbers: 98.80.Cq, 98.80.Bp, 98.80.Jk 
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' The dynamical system approach has a prominent role in modern cosmology. Any candidate to a realistic cos- 
' mological model must exhibit certain qualitative dynamical behaviors as, among others, robustness against small 
perturbations (to avoid fine tunings problems, see [l| for a recent review), specific kinds of phase space attractors 
(to describe, for instance, the recent accelerated expansion phase of the universe|3, 0] ) , and some classes of solutions 
with determined intermediate time behavior (for instance, tracker solutions |J7 candidates to explain the Cosmic 
Coincidence or Why Now? problem). In a recent workQ, Townsend and Wohlfarth present an interesting and promis- 
ing tool for the dynamical analysis of cosmological models. Essentially, they show that the equations of motion of 
fSj homogeneous and isotropic cosmological models with multiple minimally coupled scalar fields, self interacting through 
\ an arbitrary potential, do indeed correspond to the geodesic equations of a certain effective Lorentzian manifold. Such 
C ■ I correspondence, closely related to the Maupertuis-Jacobi principle of classical mechanics [g^, permits us to infer some 
' of the dynamical properties of a given cosmological model from the underlying geometry of the associated effective 
O manifold, allowing one to go a step further in the dynamical analysis of realistic cosmological models. Townsend and 
Wohlfarth, indeed, apply their own formalism to identify asymptotic accelerated expansion phases in models with 
', exponential potentials 0. As we will see, such geometrical considerations can simplify considerably the dynamical 
study of cosmological models. 

^ \ Here, we show that the work of Townsend and Wohlfarth can be largely extended, allowing the inclusion of non- 
■ minimally coupled scalar fields and also the anisotropic case, opening many new possibilities for the dynamical 
' characterization of such cosmological models. Our main results and some applications are presented in the next 



^ . sections, after the following brief introduction to the Maupertuis-Jacobi principle. 



A. The Maupertuis-Jacobi principle 



The Maupertuis-Jacobi principle in classical mechanics [Hi establishes that the dynamics of a given system can be 
viewed as geodesic motions in an associated Riemannian manifold. In order to recall it briefly, let us consider a 
classical mechanical system with N degrees of freedom described by the Lagrangian 

L{qA)^\mmi' -yiil (1) 

where i, j = 1,2, . . . ,N , the dot stands for differentiation with respect to the time t, and gij is the Riemannian metric 
on the A^-dimensional configuration space M.. All the quantities here are assumed to be smooth. The Euler-Lagrange 
equations of (jT]) can be written as 

it+^ki'i" = -9''d,V{q), (2) 
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where FJ^, is the Levi-Civita connection for the metric gij. The Hamiltonian of the system described by ([T]) 

■H{q,p)^]^g''{q)p^p,+V{q), (3) 

with Pi = gijq-' , is obviously a constant of motion, namely the total energy. For a fixed energy E, the trajectories in 
the 27V-dimensional phase-space {q^;pj) are confined to the hypersurface E = ^g^^piPj + V{q). The admissible region 
for the trajectories in the configuration space is given by 

VE = {q^M: V{q) < E}, (4) 

which can be either bounded or unbounded, connected or not. If the potential V has no critical points on the boundary 
OVe, then OVe is a iV — 1 dimensional submanifold of A4. We can easily see that if a trajectory reaches the boundary 
at a point qo, its velocity at this point vanishes and the trajectory approach or depart from qo perpendicularly to 
OVe- In particular, there is no allowed trajectory along the boundary. 

One can show that the equations of motion Q arc, in the interior ofDE, fully equivalent to the geodesic equation 
of the "effective" Riemannian geometry on Ai defined from the Jacobi metric 

g,,{q)^2{E-V{q))g.,,{q). (5) 

The geodesic equation in question is given by 

V7 d^l' , dq^ dq^ 

where u = dq^/ds is the tangent vector along the geodesic, V and F^j, are, respectively, the covariant derivative and 
the Levi-Civita connection for the Jacobi metric gij, and s is a parameter along the geodesic obeying 

^^^2{E~Viq)). (7) 

As for any classical topic, there is a vast literature on the Maupertuis-Jacobi principle. We notice only that, motivated 
by the celebrated result due to Anosov stating that the geodesic flow in a compact manifold with all sectional curvatures 
negative at every point is chaotic the Maupertuis-Jacobi principle has been recently invoked for the study of chaotic 
dynamics. (See, for instance, ^ and the references therein). 

Townsend and Wohlfarth consider homogeneous and isotropic cosmological models with N self-interacting mini- 
mally coupled scalar fields 0" taking their values in a Riemannian target space endowed with a metric Gap- The 
corresponding actions is 

'd^x^/^ (R - 5'^G„0(0)5,0"a,/ - 21/(0)) . (8) 

where R stands for the scalar curvature of the Z?-dimensional spacetime metric gij. By considering the Friedman- 
Robertson- Walker homogeneous and isotropic metric 

ds^ ^-dt^ +a^{t)dJ:l, (9) 

where I]„ represents the {D — l)-dimensional spatial sections of constant curvature n, they showed that the equations 
of motion associated to the action ([5]) do indeed correspond to the geodesies of a certain effective Jacobi (pseudo) 
metric on a Lorentzian manifold. For the spatially flat case {k — 0), for instance, the geodesies corresponding to the 
equations of motion derived from ([5]) are timelike, null, or spacelike according, respectively, ii V > 0, V — 0, or V < 0. 
Such results have been already applied to the dynamical study of the models governed by actions of the type ([S]) , see 
[H. For K 7^ 0, Townsend and Wohlfarth had to introduce some higher dimensional effective manifolds to accomplish 
their analysis. As we will see, such higher dimensional manifolds are unnecessary, all values of k can be eventually 
treated in the same framework. 

Applications of the Maupertuis-Jacobi principle to the field equations obtained from Hilbert-Einstein like actions 
have also a long history. Non- homogeneous and anisotropic cases were considered in . Applications involving distinct 
differential spaces instead of differential manifolds were discussed in [To^. Non- minimally coupled fields, however, have 
not been considered so far. 
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II. NON-MINIMALLY COUPLED HOMOGENEOUS AND ISOTROPIC COSMOLOGIES 



Non-minimally coupled scalar fields are quite common in cosmology. In particular, they have been invoked recently 
to describe dark energy (see [TTj and [l2j. and the references therein, for, respectively, models using conformal coupling 
and more general ones. See also [11, [J). The non-minimal coupling generalization of ([5]) we consider here is 

S= f d^x^ {F{cb)R ~ g''G^ii{^)d,rd,(j,f' - 2^(0)) . (10) 



Integrating by parts the action PTI)) . taking into account that R = 6H + 12H^ + Gk/o^ for the metric with 
H = a/a, we obtain the following Lagrangian 

L{a, a, (/.") = (^-m^F{^) ~ 6i/r9„F(0) + G«^((/.)0"0'5 + - 21/(0)^ (11) 

on the (A^+ l)-dimensional configuration space spanned by (a, 0"). Although we will discuss here only the D ~ A case, 
the D-dimensional one follows straightforwardly by using that R = 2{D - l)H + D{D - l)H'^ + (D - 1){D - 2)K/a^ 
in order to derive (jlip . By introducing the following Lorcntzian metric 

on the configuration space (upper case roman indices run over 1 . . . A'^ + 1), the Lagrangian (|lip can be cast in the 
form 

L{^^,^^) = GAB{<t>)^^^'' " 2Kff(</'-^), (13) 

where (^"^ = {a,(jf^) and 

V,s{4>^) = a3l/(0) - ?>KaF{4>). (14) 

It is evident the similarity between ([1]) and (fT5|) . provided that det ^ 0, which we discuss below. Before, let us 
recall that our manipulations imply that all solutions of the Euler-Lagrange of pUj) are also solutions of the Euler- 
Lagrange equations of (jlip . but not the converse. Einstein equations form a constrained system, the solutions of (jlO[) 
correspond, indeed, to a subset of the solutions of (jlll) . as one can realize by considering the Hamiltonian associated 
to ^ 

H(0^, TT^) = G^^t:a^b + 2Kff (0^), (15) 

which is a constant of motion, where it a — Gab4'^ ■ The Euler-Lagrange equations of pop . on the other hand, implies 
that 7i = (the so-called energy constraint). Hence, we must bear in mind that the relevant solutions of our original 
problem correspond, in fact, to the H — subset of the dynamics governed by pT|) . 

A proper interpretation of (|12p as a metric tensor does require the essential assumption det Gab 7^ 0. In order to 
grasp it actual meaning and implications, let us restrict ourselves to the A^ = 1 and G = 1 case, for which 



det Gab = -6a^ [F{c^) + - (F'(0))^ j . (16) 

The vanishing of the quantity between parenthesis is known for a l ong time to be associated with the existence of 
some severe and unavoidable dynamical singularities [ll| (see also [l4l]). which render the associated cosmological 
model unphysical. Here, the vanishing of ([TBI) also leads generically to an actual geometrical singularity. Thus, 
the assumption of det Gab ^ assures that the model in question is free of these singularities, geometrical and 
dynamical, a basic requirement for any realistic model. This can be viewed as the first application of the formalism; 
we will return to this issue in the last section. Note that one of the most interesting peculiarities of the conformal 
coupling [N — 1, F{(j)) — l — cjP'/Q) is that it can always evade the det Gab = singularity, since F{(j)) + 1 (F'(0))^ = 1 
for such case. 

The Lagrangian is ready to be considered under the Maupertuis-Jacobi principle. For Vcs — 0, the Euler- 
Lagrange equations of (fT3)l already correspond to geodesies of p2|) . Moreover, from the energy constraint 7i = 0, one 
has that they in fact correspond to timelike geodesies. For V^fi ^ 0, let us introduce the Jacobi (pseudo) metric 



Gab — '2-\Vcs\Gab- 



(17) 
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One can see that, in accordance with the Maupertuis- Jacobi principle, the Euler-Lagrange equations of (I13p correspond 
to geodesies of ([T7)) . parameterized by s, see Eq. ([7]). From the energy constraint, one gets 

implying that the geodesic are timelike and spacelike for, respectively, Vcs > and VeS < 0. All the results's'l that 
have motivated the present work can be obtained in a simpler way by setting F = 1. The analysis corresponding to 
the K = case has already appeared in [Tsj . 

III. THE ANISOTROPIC CASE 

Our analysis can be extended in order to include also anisotropic models. Let us consider a model with action (jlOp 
and a Bianchi I type _D-dimensional metric 

D-l 

ds^ = -dt^ + ^ aj{t)dx\ (19) 

i=l 

for which the scalar curvature is given by 

^ [h., + Hf) + H,H, j , (20) 

i—1 i—l,j>i I 

where Hi — hi/ai. Again, by integrating the action (jlO[) by parts, we obtain the fohowing Lagrangian 

L=i\{a\(~2F{^) J2 i^^i?, -2 ( ^ ff, U"a„F(0)+G„/3(</')<^>''-2n</') I (21) 

\i=l / \ i=l,]>i \i=l ) / 

on the [N + D — l)-dimensional configuration space spanned by (a^, 0"). As in the previous section, we will restrict 
ourselves here to D = 4 case, the £)-dimensional one follows straightforwardly and without surprises. By introducing 
the following metric 

(0 -a-sF -a2F -a2azdpF \ 

-a^F -aiF -aia^dpF \ , . 

~a2F ^aiF ~aia2dpF ^^^> 

-a2azdaF -aia^daF -aia2daF 0102036*0/3 / 

on the configuration space (upper case roman indices run now over 1 . . . + 3), the Lagrangian (j2ip can be cast in 
the form given by (|13p . with — (oi, 02, 03, (/)") and 

l^eff(0^) = aia2a3F(</)"). (23) 

Note that the (pseudo) metric (|22p has signature (3, A). As in the isotropic case, the determinant detG^s plays a 
central role in the geometrical and dynamical analyses. For A^ = 1 and G = 1, it reads 

det Gab = -2(010203^^(0))' (^F{cb) + ^ (F'((/)))'^ . (24) 

In this case, in addition to the singularity present for the isotropic models, we have also a singularity for F^cj)) = 0. 
Such anisotropic singularity has been already described in [14| in a dynamical analysis, where it is shown that it can 
be indeed used to rule out large classes of models. We notice, nevertheless, that the dynamical analysis of is 
considerably more involved than the geometrical analysis presented here. 

As in the isotropic case, with the introduction of the metric and the effective potential the system is 

ready to be considered under the Maupertuis-Jacobi principle. 
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IV. DISCUSSION 



The great value of the dynamical analysis comes from the possibility of ruling out large classes of models that are 
not viable from the theoretical point of view. As we have already said, any candidate to a realistic cosmological 
model must exhibit certain qualitative dynamical behaviors. For instance, the singularities corresponding to the 
condition det Gab = render the associated cosmological model inviable. Geometrically, such singularities imply 
that no geodesic can be extended beyond the singular points, what would imply a kind of future singularity in the 
associated cosmological model. This is precisely the case of the isotropic singularity of Section II and the anisotropic 
one of Section III. However, we stress that the identification of such singularities in the dynamical analysesfH, [3| 
is considerably more involved than in the present geometrical approach. For multiple scalar fields, the dynamical 
analysis is much harder (see, for instance, [l6|), and the present geometrical approach can be even more useful. 

The far most common dynamical analysis in cosmology is the stability classification of certain solutions. De Sitter 
asymptotically stable fixed points are particularly relevant to the description of the recent accelerated expansion 
phase of the universeQ. Such fixed points correspond to isotropic limits such that (j)" and H are constants, 0"^(s) = 
(006^*^"*^ 00 ). The geometrical analysis presented here can help the identification of such points. For this purpose, 
let us consider the the geodesic deviation equation, which governs the local tendency of nearby geodesies to converge 
or to diverge from each order 

0^VA</'''Vsn^ = i?^BC^7i^(/)^0^, (25) 

where Rabc^ curvature tensor of the Jacobi metric, and is a vector orthogonal to t/)'^, pointing to the 

direction of the deviation. Let us take, for instance, the same case of Section II, = 1 and G = 1, with 14ff((/)) > 
now. For the geodesic corresponding to a de Sitter fixed point, one has = {{aQH/2V^s)e"\Q). The only non 
vanishing component of (|25p reads 

n(s)=-i?i2i'^(s)(<^(^))', (26) 

The stability of the de Sitter solution requires a bounded n{s) for s — > 00 and, hence, R121' cannot be negative, 
eliminating, in this way, the possibility of finding out viable de Sitter fixed points inside regions where Ri2i^ < 0. In 
the present case, one has 



R 
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3aF 



□K=ff daV,nd^V,s FF" - \ [F'f 



\ ^"^^ [Fief.) + I {F'{<i>)Y 

For F = 1, we have -Ri2i^ — 3T^"(0o)/(2a''^F^((/)o)) on the fixed point a-nd the usual requirement of stability 
F"(0o) > is recovered. For more general cases, the criteria i?i2i^ > < '^^'^ be used to locate regions where 
de Sitter fixed points could and could not appear, respectively. We finish by noticing that our analysis is in perfect 
agreement with the comprehensive study of jl7| . 
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